Charged black holes in nonlinear massive gravity 
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We find static charged black hole solutions in nonlinear massive gravity. In the parameter space 
of two gravitational potential parameters (a, /3) we show that below the Compton wavelength the 
black hole solutions reduce to that of Reissner-Nordstrom via the Vainshtein mechanism in the weak 
field limit. In the simplest case with a — P — the solution exhibits the vDVZ discontinuity but 
ordinary General Relativity is recovered deep inside the horizon due to the existence of electric 
charge. For a ^ and /? = 0, the post-Newtonian parameter of the charged black hole evolves to 
that of General Relativity via the Vainshtein mechanism within a macroscopic distance; however, a 
logarithmic correction to the metric factor of the time coordinate is obtained. When a and P are 
both nonzero, there exist two branches of solutions depending on the positivity of /3. When (3 < 0, 
the strong coupling of the scalar graviton weakens gravity at distances smaller than the Vainshtein 
radius. However, when /3 > the metric factors exhibit only small corrections compared to the 
solutions obtained in General Relativity, and under a particular choice of /3 = /Q the standard 
Reissner-Nordstrom-de Sitter solution is recovered. 

PACS numbers: 04.50.Kd, 14.70.Kv 
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I. INTRODUCTION 



The question on whether there exits a consistent co- 
variant theory for massive gravity, where the graviton 
acquires a mass and leads to a modification of General 
Relativity, was initiated by Fierz and Pauli (FP) It 
was observed that at quadratic order the FP mass term 
is the only ghost-free term describing a gravity theory 
with five degrees of freedom [2|. However, it is not pos- 
sible to recover linearized Einstein gravity in the limit 
of vanishing graviton mass, due to the existence of the 
van Dam-Veltman-Zakharov (vDVZ) discontinuity aris- 
ing from the coupling between the longitude mode of the 
graviton and the trace of the energy momentum tensor 
[3, 01 ■ It was later observed that this troublesome mode 
could be suppressed at macroscopic length scales due to a 
nonlinear effect: the so-called Vainshtein mechanism Q. 
However, these nonlinear terms, which are responsible for 
the suppression of vDVZ discontinuity, lead inevitably to 
the existence of the Boulware-Deser (BD) ghost Q, mak- 
ing the theory unstable (7l-[l0j. 

Although for many years it was believed that the the- 
ory of massive gravity always contains BD ghosts, a fam- 
ily of its nonlinear extension was recently constructed 
by de Rham, Gabadadze and Tolley (dRGT) [H [11. 
This is a two parameter family of nonlinear generaliza- 
tion of the FP theory, where the BD ghosts are removed 
in the decoupling limit to all orders in perturbation the- 
ory through a systematic construction of a covariant non- 
linear action [l3l4l6| (see [T^l for a review). As a conse- 
quence, the theoretical and phenomenological advantages 
of the dRGT theory led to a wide investigation in the lit- 
erature. For example, cosmological implications of the 
dRGT theory are discussed in (iSmTj : black holes and 



spherically symmetric solutions were analyzed in |38l-l46j ; 
and the theory's connections to bi-metric gravity models 
were studied in HTHGOll. 



Among these phenomenological studies, is the search 
for observationally suitable spherically symmetry solu- 
tions. Theories of massive gravity can be strongly con- 
strained due to the vDVZ discontinuity appearing in the 
post-Newtonian parameters. Recently, a class of black 
hole solutions in the theory of "ghost-free" massive grav- 
ity was analyzed by Koyama, Niz and Tasinato (KNT) 
(ssl . [39| . Their result shows that the behavior of lin- 
earized solutions in General Relativity can only be repro- 
duced below the Vainshtein distance in a certain region 
of parameter space. An exact Schwarzschild-de Sitter 
(SdS) solution was found for a group of specially selected 
parameters (4l| . 

In this paper we present a family of static, electrically 
charged black hole solutions in the theory of ghost-free 
massive gravity. The solutions posses a Vainshtein ra- 
dius, below which the linearized solutions of Einstein 
gravity are approximately recovered in the weak charge 
limit. In massive gravity the longitudinal mode of gravi- 
tons is stroirgly coupled to the trace of energy momentum 
tensor and the existence of an electric charge can strongly 
affect the behavior of this mode. 

The paper is organized as follows. In Section [Hi we re- 
view the model of nonlinear massive gravity and present 
the equations of motion for gravitational fields in a spher- 
ically symmetric background. In section IIIII we investi- 
gate in detail a stellar background with a static electric 
field. In section HV] we analyze the charged black hole so- 
lutions both analytically and numerically, and we show 
that the Vainshtein effect can be made manifest in a cer- 
tain parameter space in the weak field limit. Finally, 



section IVl summarizes our results. 



II. GHOST-FREE MASSIVE GRAVITY 

Massive gravity has an effective field theory descrip- 
tion given by Einstein gravity plus the covariant FP mass 
term. For the dRGT model Lagrangian the potentially 
pathological term can be absorbed by total derivative 
terms, leading to equations of motion that are at most 
second order in time derivatives 11211. 



The dRGT action 



where arises from the graviton potential Q 
~" i i^ui, - AL/C^^ H 



+6/3 \k.% - KKl, + i/C^. { [/C]2 - [K.^] ] I 
-/?5m.{[/C]3-3[/C][/C2]+2[/C3]} . (9) 

In addition to the generalized Einstein equations, the 
Bianchi identities lead to the constraint: 



V^X„, = 
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The gravitational action is: 
1 



S 
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where R is the Ricci scalar, and U is a potential for the 
graviton which modifies the gravitational sector. Specif- 
ically, lA is given by 



(2) 



in which and are dimensionless parameters. More- 
over, IA2, U3 and are defined as 



U4 = [JC]^ - 6[/C]2[/c2] + 8[/C][/C3] 
with 



(3) 
(4) 

3[IC^f-6[IC% (5) 



(6) 



where the rectangular brackets denote the traces, namely 
[K] = /C^. Finally, in the above relation the four- form 
fields 0° are the Stiickelberg scalars introduced to restore 
general covariance Q. 



B. Generalized Einstein Equations 

For convenience we choose the unitary gauge (j)"' = 
x^(5° and thus the tensor g^^ is the observable describing 
the five degrees of freedom of the massive graviton. In 
addition, we regroup the two parameters and on of the 
graviton potential ([2|) introducing two new parameters, 
a and /?, as 



a — 1 



03 = 



«4= 2 



12 



(7) 



to simplify the background equations of motion. 

Varying the action with respect to g^^^, leads to the 
modified Einstein equations: 



(8) 



III. GENERAL ANALYSIS ON SPHERICALLY 
SYMMETRIC CHARGED BACKGROUND 

Having derived the equations of motion we now study 
the dynamics of a gravitational system described by mas- 
sive gravity under a fixed background symmetry. Follow- 
ing |38| , we consider the most general form of the metric 
respecting spherical symmetry. 



-N^{r)df + 



dr^ 



2D{r)dtdr 



r^dnl 



ill) 



where dVt\ = d9^ 



sin^ 9dip'^. 

Notice that the identity gorRoo ~ gooRor = yields 
an algebraic constraint and correspondingly leads to two 
solution branches: either D{r) ~ 0, or the effective en- 
ergy density is proportional to gu- The latter case always 
yields a constant metric factor H{r) with value depend- 
ing on specific backgrounds (4^ . In the following we will 
focus on the first case D{r) = 0, corresponding to a di- 
agonal metric. 



A. Charged spherical symmetric background in 
General Relativity 

We consider a generic Maxwell field in curved 
spacetime, with standard Lagrangian. The Maxwell 
equations are 



-gFn 



(12) 



where J'^ is the current density. For a static electric 
charge Q in the gravitational system, the components of 
the Maxwell field are: 



For = E{r) , Ee 



E,, 



, B = . (13) 



^ Note that the Einstein equation shown in |44|| contains a typo of 
an extra 1/2 factor, however the equations of motion which give 
rise to solutions were based on the method of varying the action 
with respect to the metric factor directly. 
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With vanishing source term, the inhomogeneous 
Maxwell law gives 



yielding the solution: 

QNH 







Eir) 



(14) 



(15) 



where Q is an integration constant which is typically in- 
terpreted as the electric charge. The factors N, F, and H 
are the fields introduced in the background metric ([TT|). 
Varying the electromagnetic Lagrangian with respect to 
the metric gives the energy momentum tensor: 



(16) 



It is well known that a static and spherically symmetric 
solution under standard General Relativity is described 
by the Reissner-Nordstrom (RN) solution: 



with 



- 1 



r^dn\l7) 



An 



rs = 2GM 



(18) 



being respectively the length scale associated with the 
electric charge Q and the Schwarzschild radius deter- 
mined by the mass of the spherical object M. The above 
solution exhibits a singularity at r = 0, in which the met- 
ric coefficient becomes zero and invariants like the Ricci 
and Kretschmann scalars diverge, however it is usually 

shielded by a horizon at r = (^rs + y^r| — 4rQ^ /2, in 

which the metric coefficients become zero in the above 
coordinates but the invariants remain finite, thus satisfy- 
ing the cosmic censorship and no-hair conjectures. Note 
however that negative mass or highly charged solutions 
do not have a horizon, and thus a naked singularity ap- 
pears. 



B. Equations of motion for splierically symmetric 
charged background 

Going beyond the above General Relativity solution, 
in the case of nonlinear massive gravity, we must con- 
sider the effects of the graviton potential. Thus one can 
combine the generalized Einstein equations ^ , the solu- 
tion to the Maxwell field (|15p , and the energy momentum 
tensor p6p , and obtain three main equations of motion. 

We consider the spherically symmetric metric Ansatz 
PTjl with D{r) =0. The background equations of motion 
are derived from the generalized Einstein equations ([5]). 



After some algebra, we extract the equations of motion: 
For the "00" component of generalized Einstein equation: 



(1 + m^r^)H^ + 2m^r^{F - 3)H^ 



GQ^H^ 

-H^[2r{FF ~ SmV^) + 3m^r^F + F'^] 
+2rFH[F{rH + 3ij) + rFH\ - br^F'^H^ 
+m^r^{H - 1)H^[2H{1 - a + 3^) - 6/3 
+F{1 -a + 6(5- 3H{1 - a + 2/3))] . (19) 

The "rr" component of generalized Einstein equation 
takes the form: 

^^'^5^ = (1 + m\^)NH^ + 2m\\l - 'SN)H^ 
47rr^ 

+H'^[N{6m^r'^ ~ F^) - ri2F^N + 3mV^)] 

+2rF^HH{rN + N) - r^F^NH^ 

+m^r^{H - 1)H^[1 -a + 6/3(1 - N) 

-H{3{1 - a) - 2N{1 - a + 3/3) + 6^)] . (20) 

Further, the generalized Einstein equation along the solid 
angle is given by. 



= m^rH^[{3 - F)N - 1] + 2rF^NH^ 



GQ^NH^ 

-FH[rNFH + F{rNH + rNH + 2NH)] 
+H'^[rFFN + NFF + F^{rN + N) 
+m^rF{3N - 1) - Gm'^rN + 3mV] 
+m^rH^{{l - a)[4:N - 3 + H(2 - 37V) 
+F{2-SN + H{2N - 1))] 

+2{F -1){H- \){N - 1)(2 ~2a + 3/3)} . (21) 

The equations of motion may also be obtained by vary- 
ing the action with respect to the metric fields N . F and 
H , respectively. In addition, there is another constraint 
equation from the Bianchi identity ([TU]): 
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f{ H[2 - "irN - 2iV(3 + rH)] 



rNH 

-2(3iV - l)rH + 2H^{rN + A^) I - 2H^[l + N{H - 3)] 



+ (1 - a)|2i/2[2 - 3A + H{2N - 1)] 
+F[rNH^ + 2rij(2 - 3A^) + H'^{2 -AN ~ ArN) 
+H{6N -A- 2rH + 3rN + ArHN)]^ 

+2{2a - 3/3 - 2){H - 1)I^2H^{N - 1) 

+F[rNH^ + 2rH{N - 1) - H{2N -2 + rN)] 1 1 , (22) 



where in the above equations the dot denotes a derivative 
with respect to the radial coordinate r. 

Note that when we take Q — the above equations of 
motion reduce to the case of the strong interactive solar 
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system governed by nonlinear massive gravity, discussed 
in (39| . Our expressions are in agreement with theirs, 
expect that the radial metric factor is in our case but 
becomes F in their convention. The relations between the 
parameter spaces (03,04) and {a,/3) are already given 
in ([7]). Equation ([22|) obtained from the Bianchi identity 
constraint is not an independent equation after we choose 
D{r) = 0. 



nonlinear massive gravity might be viable only when the 
invariant is non-singular. Though its divergence does 
not bring any manifest problem to observable variables 
at background level, it may be a problem for fluctuations 
passing through the horizon. In the present paper, we 
focus on the dynamics of background solutions outside 
the horizon and thus we do not address this issue. 



C. The linearized treatment of nonlinear massive 
gravity 

Following the idea developed by KNT [H [Hi, we 
study the solutions to such a gravitational system in 
the weak field limit. We can expand the metric factors 
around a Minkowski background as 



N{r) = 1 + n{r) , 
F{r) = 1 + f{r) , 
H{r) = 1 + h{r) , 



(23) 



and then investigate the linear perturbations. However, 
we need to be aware of the fact that the factors n and / 
can be treated as linear perturbations as in General Rel- 
ativity, while h could, in principle, take large values since 
this factor corresponds to the strong interactive nature of 
the scalar mode of graviton in solar system. Therefore, 
we need to keep higher orders in h and truncate equations 
of motion to leading order of n and /. We demonstrate 
this behavior both in analytical and numerical calcula- 
tions in the following. 

Before expanding the background equations perturba- 
tively, we rescale the radial coordinate by introducing a 
new metric variable 



r 

H 



and correspondingly introduce a new metric factor 

1 + ./ 



1 + / = 



l + h + ph' 



(24) 



(25) 



where the prime denotes a derivative with respect to p. 
As a consequence, the linearized metric can be expressed 



-[1 + 2ji{p)]dt^ + [1 - 2/(p)]<ip2 + p2^j^2 (-2g) 



which is asymptotic to Minkowski background when n 
and / become negligible. 

Apart from the usual curvature invariants, nonlinear 
massive gravity presents a new basic invariant incor- 
porating both the metric and the Stiickelberg scalars, 
namely 1°'^ = g'^'^dfj,(j)°'d^(j>^ . Under unitary gauge it 
seems that this invariant encounters a divergence on the 
event horizon if one takes the Minkowskian asymptotic 
at large scales (6^ . that is, one obtains singularities in 
the place where General Relativity had simple horizons. 
The authors of ji^l argue that a black hole solution to 



D. The post-Newtonian parameter and the vDVZ 
discontinuity 

To study the nontrivial effects of massive gravity we 
consider the metric factor h as a perturbation mode. By 
expanding the equations of motion ([T9)) . (|20|) and ((22|) 
presented in the Appendix, consistent with the linearized 
background (|26p . we get 



(2 + mV^)/ + 2p/' + m^p^{3h + ph') 
V^(2/i-n) = 
pn' + 2/ = , 



2pn' + 2/ + ?7iV'(2/i -n) = -^{l + n) 



,(27) 
(28) 
(29) 



up to leading order. 

Under the above approximation the metric factors are: 



n{p) 



'Sp 
GmQ^ 
IGnp 



[e™^Ei(-77ip) - e-™''Ei(mp)] , (30) 



fip) 



2GMe-'^P{l + mp) GQ^ 



3p 



87rp2 



X [(1 - mp)e"^PEi{-mp) 
-{1 + mp)e-"'PF.i{mp)] , 



(31) 



and 



h^p)-^ ———{l + mp + m^p^) + 



+ 



SmPp'^ 
GQ^ 
32mnp^ 

-(l + mp + m2p2)e-™PEi(mp)] , (32) 



167rp2 

X f(l - mp + m^p^)e"'PEi{-mp) 



where "Ei" is the exponential integral function defined 
by Ei(x) = J'^^e* dint. 

Let us now examine the post-Newtonian parameters 
for nonlinear massive gravity. The post-Newtonian pa- 
rameters are strongly constrained by solar system obser- 
vations (see (6l| for a detailed introduction and a com- 
prehensive review), and therefore may be used to place 
strong constraints on a given theory. The first post- 
Newtonian parameter 7 is defined as the ratio of / and 



/ 

n 



(33) 
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in the weak field limit. In the case of General Relativity 
7 = 1, as we can immediately read from the RN solu- 
tion (fT7|) . However, in the regime of length scales much 
smaller than the Compton wavelength (/9,„ = I/'ti) in 
nonlinear massive gravity, this parameter can be approx- 
imately expressed as 



7 ■ 



327rA/p(l + mp) - 6Q^e' 
647rMp- GQ^gmp 



(34) 



where we have made the reasonable assumption that the 
solar-system mass M is much larger than graviton mass. 

It is straightforward to see that when the solar sys- 
tem does not carry an electric charge ((5 = 0), we obtain 
7 ~ (1 + mp)/2, and thus in the massless limit it reduces 
to 7 = 1/2, which is in stark disagreement with the value 
in General Relativity. This behavior exactly manifests 
the vDVZ discontinuity, and thus we face the difficulty 
of explaining solar-system observations at present. How- 
ever, when the contribution of Q is taken into account, 
the post-Newtonian parameter could approach to 1 when 
p is much smaller than the length scale 



Pg 



32nM 



(35) 



This result implies that the vDVZ discontinuity exists 
in a charged solar system governed by nonlinear mas- 
sive gravity, but the explicit behavior is different from a 
electro-neutral one. 

Finally, from the solutions ^ and ^ we find 
that the weak-field limit is a good approximation at large 
distances. However, when the radial coordinate p de- 
creases below a critical value 



Pv 



fGM\ 



1/3 



(36) 



the absolute value of the factor h will increase exponen- 
tially and become much larger than unity. This critical 
radius is the so-called Vainshtein radius. We perform a 
detailed analysis on the perturbation equations by keep- 
ing all nonlinear order terms of h in the next section. 



IV. CHARGED BLACK HOLES AND 
VAINSHTEIN MECHANISM 

Due to the famous Vainshtein mechanism the scalar 
degree of freedom in massive gravity becomes strongly 
coupled in the limit of small graviton mass, and thus the 
linearized treatment performed in equations (P7|) . (pS)) 
and is no longer reliable. This is observed by follow- 
ing the evolution of metric factor h below the Vainshtein 
radius, where the absolute value of h becomes much 
greater than unity. Consequently, although we can keep 
treating n and / as small perturbations in this regime, 
higher order terms in h should be taken into account. 



A. Perturbation equations with nonlinear 
corrections 

Keeping leading order in n and /, the perturbed equa- 
tions of motion including nonlinear corrections of h are 
given by 



2/ + 2p/' 



2 2 

m p 



{1 - 2ah + 6l3h^)[{l + f)ph' 



+{l + h)f] + 3h{l -ah + 2l3h^) 



GQ' 
47rp2 



2/ + 2pn' - m^p^ 



n — 2(1 + n + an)h 



+{a - 



Q(in)h^ 



GQ2(1 



, 



0, (37) 

pn [I - 2ah + eph^] + 2f[l - ah] ^ , (39) 

of which the first two correspond to the (00) and (rr) 
components of the generalized Einstein equations and the 
last one arises from the perturbed Bianchi constraint. 

We first solve equation ([37|) by neglecting all high-order 
terms proportional to m^, G^M^, m?GM, G^Q"^, and 
m?GQ^ . Therefore, the metric factor / can be expressed 
in terms of h: 

GQ^ GM 
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m^p^ 



87rp2 



{h - ah^ + 2/3/i^) . (40) 



Inserting the expression (|40p into the perturbation equa- 
tion (|38p . we obtain the radial derivative of n as a func- 
tion of h: 

GQ^ GM 



irP'p 



Anp^ 



{h - 2/3/1'^ 



(41) 



and thus the metric factor n can be acquired by perform- 
ing integration. The key to solving for the metric factors 
71 and / is to extract the solution for h. Therefore, we 
combine the expressions PO)) . (HI]), and the perturbed 
Bianchi constraint p9p. and then we derive the polyno- 
mial equation: 



2 2 
TO p 



-h 



4:TTp 



3ah^ 



(a^ +4/3) /i^ + 6/3^/1^ 



(42) 



in which all nonlinear terms of h have been taken into 
account. 

After having obtained equations (gO]), (|4T]) and (|42| . 
we are able to calculate the explicit forms of the metric 
factors n, f and h under different parameter choices of 
(a,/3). In the following subsections, we investigate these 
solutions further. 



B. Case I: a = /3 = 

Let us first consider the special subclass with a = /3 = 
0. In this case all higher order terms of h vanish automat- 
ically, and thus our task reduces to solving the linearized 
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equations (|27p . ([^5]) and ([^U)) . Therefore, the correspond- 
ing sohitions of metric factors are aheady given in equa- 
tions ([30 | . ([3T|) and (|32l) . 

The post-Newtonian parameter takes the value of 1/2 
in the regime of length scales of interest. Although this 
parameter increases to 1 at very small values of radial 
coordinate due to the effect of electric charge, the cor- 
responding length scale is deeply inside the horizon of 
the black hole. As a consequence, the standard General 
Relativity result can not be recovered in the solar sys- 
tem governed by massive gravity with such a parameter 
choice, and therefore this case is already observationally 
ruled out. 



C. Case II: Q / and /3 = 

We now proceed to the case with a vanishing /3 but non- 
vanishing a. In principle, one can solve the /i-equation 
P9p exactly, however the resulting expression is quite 
complicated, hiding the underlying physics. Therefore, 
instead of finding an exact solution to h we solve the per- 
turbed Bianchi equation approximately by keeping domi- 
nant terms in h. Different from the electro-neutral case of 
[3^ , both the solar-system mass M and its electric charge 
Q contribute to the L.H.S of equation ()42|) . Therefore, 
we need to solve (j42p by assuming its L.H.S is dominated 
by M and Q respectively. 

We first consider that the contribution of the electric 
charge Q becomes dominant in the L.H.S of (j42|) when 
p <ti pv- In this limit we solve equation (|42p and expand 
h up to order 0{p^): 

h - 



3 

2^ 



PyP 



(43) 



2^Pq ' 

where we have introduced a new parameter for a critical 
length scale 

1/4 



PQ 



\47rm2/ 



(44) 



One can see that the approximate expression for h is valid 
only when p < pg. Comparing with the non-charged 
system analyzed in (39| . we find that the leading term 
of /i in a charged solar system is proportional to p^^ 
instead of p^^. Moreover, the suppression scale is pg 
rather than the Vainshtein radius pv- However, since 
for a group of canonical parameters accommodating with 
astronomical data pv is usually much bigger than pg, 
we deduce that the contribution of the third term in the 
R.H.S. of equation P5)l is considerable in a wide regime 
of length scales. Then, from equations (PU]) and (PT|) we 
can derive the expressions for n and / through a Taylor 
expansion 



GQ2 


GM 
h 


m^Pq 


87rp2 


P 


2ai/2 


GQ2 


GM 


m^Pq 


47r/92 


P 


Q,l/2 



In(mp) , 

GMp 
^ 2ai/2p2 



(45) 



(46) 



up to order 0{p'^). By observing and (|46p . when 
2GM < p < pq, the corrections of massive gravity to 
the General Relativity results are quite small. When p 
is larger than pq but smaller than pv, the higher order 
corrections are suppressed by py and the formulae are in 
agreement with Ref. (39| . 

When the radial coordinate p evolves to the regime 
which is close to pq but still smaller than py, the main 
contribution of the L.H.S of (|42p comes from the mass 
term M . In this case, we only keep the leading order in 
h, obtaining 



Pv 





Q;2/3p 
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n ~ 


GQ2 


GM 


GMp 


87rp2 


P 


^ 2a2/3p^ 




Gg2 


GM 


GM 


87rp2 


P 


^ 2a^/^pv 



GMp 
2a2/3p2, 



(47) 
(48) 
(49) 



Comparing our results with the analysis of [39| we find 
that the above three solutions are consistent with those 
in the solar system without a charge. This implies that 
there must exist an intermediate regime along the radial 
coordinate in which the behaviors of the metric factors 
of charged black holes are the same as those of neutral 
black holes. 

In order to provide a more transparent picture of the 
dynamics of the charged solar system described by non- 
linear massive gravity with the above parameter choice 
(a 7^ and /3 = 0), we numerically evolve the perturba- 
tion equations. In particular, without loss of generality 
we consider a = 1, and we choose the graviton mass 
m = 10~20Mp and the solar-system mass M = lO^Mp, 
setting also the Planck unit Mp = 1/G = 1. Correspond- 
ingly, this group of parameters yield the Compton wave- 
length pm ~ 102'', which is much larger than the Vain- 
shtein radius pv = 2.15 x 10^^. Moreover, we consider 
the source term of the matter field to be a weak charge 
Qw = 1-77 X 10^ or a strong charge Qs = 3.54 x 10'', 
respectively. 

In Fig [T] we depict the absolute value of the metric 
factor /i as a function of the radial coordinate p. From 
the above parameter choices we find that the suppression 
scales associated with the electric charges are given by 
Pqw = 2.24 X 10^^ (which is represented as a purple 
dashed line for the weak charge Qw) and pqs = 10^'^ 
(which is represented as a pink dashed line for the strong 
charge Qs)- The wide sparse shadow regime denotes the 
space between the inner horizon and outer horizon in the 
case of the weakly charged solar system, and the narrow 
dark shadow regime denotes the space between the two 
horizons in the case of the strongly charged solar system. 

From Fig. [l]we that there exist three different slopes 
for the metric factor h along the radial coordinate p in a 
charged solar system. When p is greater than the Vain- 
shtein radius pv but less than the Compton wavelength 
Pm, \h\ scales approximately as p~^ and its value is much 
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Figure 1. Plot of the evolutions of the absolute value of metric 
factor h as functions of radial coordinate p in a charged so- 
lar system described by nonlinear massive gravity. The model 
parameters are taken as: a — 1 and (3 = 0. In the numerical 
calculation, we take m = 10^^" and M = 10^. The corre- 
sponding Compton wavelength pm = 10^° and the Vainshtein 
radius pv = 2.15 X lO'^^ are denoted on the top of the figure. 
The values of weak and strong electric charges are provided m 
the plot. All dimensional parameters are in Planck units. 



smaller than 1 , and therefore the weak field limit at large 
distances is valid. When p evolves to be smaller than pv, 
\h\ becomes much larger than unity quickly and scales 
as p~^ . This behavior is in precise agreement with (|47p . 
When p decreases to a regime which is much shorter than 
pg, we observe that the slope oi\h\ changes again which 
gives rise to \h\ ~ . This transition on |/i|-slope im- 
plies that the contribution of electric charge Q becomes 
dominant in determining the dynamics of scalar graviton. 
Moreover, the length scale for the slope transition on h 
depends on the value of pQ and thus it is determined by 
the combination of Q and m. For a fixed graviton mass 
TO, the value of pq in a strongly charged solar system (de- 
noted by the purple dashed line) is much larger than that 
in the case of weak charge (denoted by the pink dashed 
line). Correspondingly, the transition of |/i|-slope in the 
case of strong charge occurs at a larger distance as shown 
in the red solid curve, while the transition of |/i|-slope in 
the case of weak charge takes place at a smaller distance 
as shown in the blue dashed curve. 

In Fig. [2] we depict the ratios n' /n'cR and f / fcR as 
functions of radial coordinate p. From the upper graph 
we observe that the correction to the metric factor n from 
the electric charge Q is very small, since the three curves 
(a red solid line representing for strongly charged case, a 
blue dashed line representing the weakly charged case, 
and a black dotted line representing the non-charged 



case) almost coincide. Additionally, when p is outside 
the Vainshtein radius, the ratio of n' in nonlinear mas- 
sive gravity and n' in General Relativity takes the value 
of 4/3. This numerical result is in agreement with the 
analytic result obtained in the previous section when one 
compares the first term in the R.H.S of equation (|30p and 
the second term in the R.H.S of equation (|^5|) . However, 
this ratio tends to 1 when p lies in the regime oi p < pv 
due to the Vainshtein effect. As a consequence, the met- 
ric factor n of nonlinear massive gravity roughly recovers 
the General Relativity result inside the Vainshtein radius. 



PQW Pqs Pv P, 




10"' 10' 10' 10' 10' 10' 10" 10" 10" 10" 10" 
P 



Figure 2. Plot of the evolutions of the ratios n /n' qr and 
f / fcR o,s functions of radial coordinate p in a charged solar 
system described by nonlinear massive gravity. In the numer- 
ical calculation, the parameters are chosen to be the same as 
those provided in Fig. QJ The "f" m the lower panel repre- 
sents for the metric factor f in the main text. 



From the evolutions of / in charged gravitational back- 
grounds, depicted in the lower graph of Fig. [21 we can see 
that they are different from the non-charged one when p 
approaches to the inner horizons. This yields the non- 
trivial modification to the post-Newtonian parameter at 
small length scales as shown in Fig. |31 When the radial 
coordinate is smaller than the inner horizon of a charged 
black hole we get 7 = 2, which disagrees with the Gen- 
eral Relativity case 0. From Fig. El one can see that 
inside the Vainshtein radius (but outside the outer hori- 
zon) we obtain approximately 7=1 and thus such a 



^ However, one needs to be aware of the fact that in this regime the 
perturbative treatment of n and / is not valid and thus a com- 
pletely non-perturbative analysis is required. Such an analysis 
lies beyond the scope of the present work. 
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solar system could conform with observations. When p 
becomes greater than py, the post-Newtonian parameter 
7 evolves to 1/2 which is expected by the effect of vDVZ 
discontinuity as shown in p4p . 




10' 10' 10' 10' 10' 10' 10" 10" 10" 10" 10" 
P 



Figure 3. Plot of the evolutions of the post-Newtonian param- 
eter 7 as functions of radial coordinate p in a charged solar 
system described by nonlinear massive gravity. In numerical 
calculation, the parameters are chosen to be the same as those 
provided in Fig. QJ 



1. /3 < 

Let us first consider the subcase j3 < 0. By solving 
h numerically we find that there is only one branch of 
solution which is real between the Schwarzschild radius 
and infinity. It evolves from a very large value to a small 
constant, as moving from the Schwarzschild radius to the 
Compton wavelength. Around the Vainshtein radius it is 
of order 0(1). A strong charge suppresses the value of h 
slightly when it is close to and inside the Schwarzschild 
radius. Beyond this range the charge does not cause an 
obvious effect. This can be seen in Fig. SI where the 
three curves almost coincide. 



Q=1. 77x10' 

Q=3.54x10° 

■ ■ ■ ■ without Q 
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D. Case III: /? / 

The case where both a, (3 parameters are non vanish- 
ing can be divided in two subcases: /? < and (3 > 0. 
Generally, we are unable to obtain the exact solution of 
the metric factor h from equation (f^^ . Therefore, we 
can only solve the equations of motion semi-analytically 
by keeping the leading order terms and then compare 
with numerical computations. However, it is interest- 
ing to notice that, in the case of /3 > there exists 
a class of exactly analytic solutions for a special fam- 
ily of parameter choice /? = which was also ob- 
served in Schwarzschild-like solutions in massive gravity 
mil [Hi [13 1- In the following we will analyze these cases 
in both analytical and numerical ways, respectively. 



^ An RN solution was constructed in the de Sitter (dS) background 
in [4ll . |44|| . In the present work we solve the exact equations 
of motion of a charged solar system in massive gravity and we 
obtain an exact form of RN-dS which will be shown in later 
context. 



Figure 4. Plot of the evolutions of metric factor h as functions 
of radial coordinate p in a charged solar system described by 
nonlinear massive gravity. The parameters of the massive 
gravity model are taken as: a — 1 and jS = —1/2. Moreover, 
m and M are the same as those provided in Fig\^ 



As shown in this figure /i >> 1 in the range rs < p < 
Pq. Thus, in the semi-analytical calculation, we keep the 
leading order terms of the equations of motion for the 
metric factors, and then we acquire 



h : 



(Pol _ 



3U/3 



2/31/3^ p 

GQ^ am^ pq'^ 3 2/3 



2/31/3 



(50) 
(51) 

(52) 



From the above semi-analytic results we find that the 
corrections to / and n are so dramatic that the usual 
Schwarzschild-like gravitational potential (in form of 1/p) 
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is exactly canceled. This result agrees with the conclu- 
sion of [39| in which a neutral solar system was consid- 
ered. 

Finally, for the gravitational potential associated with 
the electric charge, the RN-like factor ^ in gu also disap- 
pears, but the dominant term in the small radius regime 
roughly takes the form of p^^^. For the grr component 
the sign in front of ^ changes from positive to negative 
compared to the General Relativity result. These new 
features suggest that the charged solar system described 
by nonlinear massive gravity is much more smooth near 
the origin than that described by General Relativity. 
However, since below the Vainshtein radius the differ- 
ence from General Relativity is quite significant in this 
case, the corresponding parameter space is likely ruled 
out by local, solar system experiments. 

2. /? > 

We now consider the subcase (3 > 0. Apart from the 
previous solution in which h is large-valued below the 
Vainshtein radius, there exists a second solution in which 
the metric factor h always takes a small value outside the 
Schwarzschild radius. Similarly to the previous subcase 
we solve the equations of motion numerically, obtaining 
the solution for the metric factor h as shown in Fig. [51 In 



that of General Relativity, in Fig. [S] we plot the ratios 
"V'^Gi? a'^d f/fcR, respectively. 
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Figure 6. Plot of the evolutions of the ratios n'/nQn, f/foR, 
and the quotient j = ^ as functions of radial coordinate p in 
a charged solar system described by nonlinear massive gravity. 
The parameters of the massive gravity model are taken as: 
a = 1, P — 3. Other parameters are the same as those used 
in Fig. The 'f ' in the plot represents the metric factor f 
in the main text. 




10' 10' 10' 10" 10" 10" 10" 10" 10" 

P 



Figure 5. Plot of the evolutions of metric factor h as func- 
tions of radial coordinate p in a charged solar system described 
by nonlinear massive gravity. The parameters of the massive 
gravity model are taken as: a — 1 and P — 3. Other parame- 
ters are the same as those used in Fig. 

addition, we extract the evolutions of / and n' along the 
radial coordinate p. In order to provide a clearer picture 
of the difference of the dynamics in massive gravity from 



From Fig. [5] we can see that in a strongly charged solar 
system the absolute value \h\ can become larger than 1 
inside the Schwarzschild radius. As we move away from 
the Schwarzschild radius h evolves to a constant, which 
coincides with the value when the charge is weak. After 
crossing the Vainshtein radius h approaches rapidly. By 
comparing / and n' with the results of General Relativity, 
we can clearly see the Vainshtein mechanism from Fig [51 

Since \h\ is always smaller than unity outside the 
Schwarzschild radius, we can solve its equation of mo- 
tion by neglecting the terms proportional to and h^. 
Correspondingly, the approximate solution is given by 

12(ap4 _ 2f3pQi + 2pppv^) 

\ _ r A%ppv\ap^ -2Ppq'' + 2pppv^) ^"^\ 
\ [ i3p^-2apQ^y \ j ' 

(53) 

which proves to be in good agreement with the numer- 
ical solution. The corresponding expressions for / and 
n' are obtained from equations (UHl) and (HI]), however 
they are quite complicated and thus we do not present 
them explicitly. Furthermore, we would like to point out 
that in this case the dynamical features of the metric fac- 
tors / and n are quite similar to the case of a 7^ and 
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/3 = 0, as can be seen by comparing Fig. [S] and Fig. (3] 
This implies that the solution obtained in the (3 = case 
might be dynamically stable in the parameter space of 
the nonlinear massive gravity model. 



3. 13 = o? 1%: An exact analytic solution 

In last subsection we examine a family of solutions in 
the model of nonlinear massive gravity under the partic- 



ular parameter choice /3 



/6. This special parameter 



choice was first noticed in |41l |. where the exact SdS and 
RN-dS solutions were obtained with an arbitrary cosmo- 
logical constant term. This relation was also applied by 
the authors of [i^l, who constructed a special family of 
black hole solutions on a fixed dS background. 

We insert the relation (3 = into the nonlinear 

equation of motion and we find that there exists a 
special solution for h, namely 



a 



(54) 



which implies a constant metric factor H ~ (1 + a) /a. 
Then we substitute the solution (|54p into the exact back- 
ground equation of motion ([5]) and the Bianchi constraint 
([22| . It is easy to verify that the constraint equation ((22|) 
is automatically satisfied. Working with the r coordinate 
directly we see that the main equation ([5]) yields the fol- 
lowing exact solution: 



(l + a)2 GQ^il + aY 



- ^ - (55) 
r 6a 



where tm is an integration constant. 

We can perform the following coordinate rescaling: 



t 



-t 



1 



l + a 

and introduce two coefficients 



{1 + aY 



'3a 



(56) 



(57) 



which are related to the Schwarzschild radius and the de 
Sitter radius, respectively. We get the exact form of the 
RN-dS-like solution as 



-A{r)dt^ 



Air) 



r^dn"^ 



with 



A{r) = 1 + 



r 



(58) 



(59) 



where the forms of rg, rs ^'iid r\ are provided in ([T^ 
and ((57|) . respectively. The above solution can recover 
the standard RN result in General Relativity, and fg co- 
incides to the usual Schwarzschild radius rs when we take 
m = 0. Furthermore, our result is in agreement with the 



one obtained in [4l|, however we did not introduce an 
additional cosmological constant in order to see whether 
and how a pure massive gravity model can yield a dS 
background by itself. 

Moreover, the coefficient r\ is determined by the com- 
bination of the graviton mass m and the model parameter 
a. When a is of order 0(1) we can apparently observe 
that the property of dS background is completely deter- 
mined by the graviton mass. If we apply this feature 
in a cosmological setup, the existence of a tiny graviton 
mass can drive a late-time acceleration of the universe 
and thus might explain the present cosmological observa- 
tions. However, if we tune a to an extremely large value 
then the effect of the graviton mass can be decreased 
and the background dynamics approach those of General 
Relativity. Mathematically this effect can be seen by the 
way /C appearing in the graviton potential U is roughly 
proportional to 1/a and thus it yields the amplitude of 
the effective energy-momentum tensor X^^ in form of 
1/a. Eventually, the behavior of nonlinear massive grav- 
ity approaches standard Einstein gravity in the limit of 
a > 1. 



V. CONCLUSIONS 

In the present work, we investigated the spherically 
symmetric solutions of a charged solar system in the con- 
text of nonlinear massive gravity. Due to properties of 
the graviton potential in the dRGT model, the BD ghost 
which historically plagues all massive gravity theories can 
be removed. However, inherited from other massive grav- 
ity models, the longitudinal mode of gravitons is strongly 
coupled and thus greatly affect the gravitational potential 
at macroscopic scales. Therefore, this model is expected 
to be constrained by solar system observations. 

Depending on the different dynamics of our solutions, 
the solution parameter space can be roughly categorized 
into three parts: 

• The first class corresponds to a = /3 = and thus 
the graviton's potential takes a fixed form. The so- 
lution in this subclass is well described at the per- 
turbative level, but the vDVZ discontinuity cannot 
be avoided. However, the post-Newtonian parame- 
ter in this class fails to agree with General Relativ- 
ity and thus the corresponding parameter choice is 
observationally ruled out. 

• In the second subclass, we keep /3 = but we allow 
a to be an arbitrary constant. The corresponding 
solution shows that General Relativity can be re- 
covered between the outer horizon of the black hole 
and the Vainshtein radius by virtue of the Vain- 
shtein mechanism. This scenario is similar to the 
case of the neutral black hole in massive gravity. 
However, the existence of an electric charge could 
increase the value of the metric factor h within 
a newly defined radius pq and thus the detailed 
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evolutions of time-like and space-like metric com- 
ponents behave differently from those of a neutral 
black hole. Namely, the metric factor n obtains a 
logarithmic correction when the radius is close to 
the outer horizon. 

• The third subclass of parameter choice is the most 
general in the parameter space, which requires both 
a and /3 to be non- vanishing. In this case the dy- 
namics of solutions behave dramatically different 
depending on the positivity of /3. When /3 < the 
strongly coupled scalar graviton greatly decreases 
the strength of gravity at small length scales, and 
thus the usual Schwarzschild-like gravitational po- 
tential totally disappears which severely challenges 
all astronomical observations. However, if (3 is 
positive General Relativity can be recovered again 
through the Vainshtein mechanism. This behav- 
ior is similar to the solution in the second subclass 
with /3 = 0. Therefore, the solution in this case, 
together with the solution in the second subclass, 
might provide a certain parameter space for nonlin- 
ear massive gravity to conform with current solar 
system observations. 

• Finally, there exists a particular parameter choice 
in the last subclass which suggests /3 = Un- 
der this condition the background equations of mo- 
tion can be solved exactly and yield a solution 
which is identical to the RN-dS form in which only 
the dS radius r\ contains the model parameter a. 
The exact solution with such a special parameter 
choice can recover the standard result in General 
Relativity in the limit of either a vanishing graviton 
mass or an extremely large value of the parameter 
a. 

Note that, the structure of solutions does not only de- 
pend on the sign of /3, but also on a diverse structure in 
the parameter space of a and (3. The authors of [i^l ana- 
lyzed the whole parameter space of a neutral solar system 
under nonlinear massive gravity and found that each re- 
gion showed a completely different behavior with respect 
to the inner solutions (near the body) and the asymptotic 
solutions, and thus affect the Vainshtein mechanism. It 
would be interesting to perform a global analysis on the 
parameter space in our case too. 

Lastly, we would like to mention that in the present 
work we only focus on the analysis of background so- 
lutions of a charged solar system in the dRGT model, 
without studying the perturbations. At the background 
level the parameter space of model parameter has already 
shown rich behaviors and a sizable regime has already 
been ruled out by observations. Due to the Vainshtein 
mechanism there exists an island in the parameter space 
which is consistent with astronomical data at the back- 
ground level. Moreover, this situation could dramatically 
changed if the perturbations were taken into account and 
we expect the model parameters appearing in the nonlin- 



ear massive gravity would be further constrained. How- 
ever, we leave such a study for future investigation. 
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